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. Introduction  and  Overview 


In  practice,  long  life  items  are  often  subjected  to  larger  than  normal 
stresses  (doses)  in  order  to  obtain  failure  data  in  a short  amount  of  test 
time.  It  is  common  to  test  the  items  at  more  than  one  stress  level,  and 
I typically,  more  items  are  tested  at  the  higher  stress  levels  than  those  at 

the  lower  stress  levels;  such  tests  are  known  as  accelerated  ot  overetreee 
life  tests.  The  basic  aim  is  to  make  inferences  about  the  life  distribution 

.Vi'  1 

of  the  items  at  the  normal  stress  levels  using  the  failure  data  from  accelerated 

. J 

tests. 

The  current  approach  to  this  important  problem  makes  inferences  under 
parametric  assumptions.  This  may  be  valid  in  some  situations,  and  does  yield 
results  which  are  appealing  from  a statistical  point  of  view.  However,  an 
engineer  or  a statistician  working  under  less  well-defined  conditions  may  find 
this  approach  too  restrictive.  For  example,  it  is  common  to  assume  that  at  all 
stress  levels,  failure  times  are  governed  by  exponential  distributions  or  Weibull 
distributions.  In  addition,  a functional  relationship  between  the  parameters 
of  the  failure  distribution  and  the  applied  stress  is  assumed.  Such  a relation- 
ship is  known  as  an  acceleration  function  or  a time  transformation  function ; 
examples  can  be  found  in  Mann,  Schafer,  and  Singpurwalla  (Chapter  9,  (1974)). 

In  two  recent  papers,  Phatak,  Zimmer,  and  Williams  (1977)  and  Shaked,  Zimmer 
and  Ball  (1977),  the  distributional  assumption  is  dropped,  but  the  requirement 
is  retained  that  the  (unknown)  failure  distribution  be  of  the  same  form  at 
the  use  and  at  all  the  accelerated  stress  levels. 


1 


Clearly,  there  are  many  situations  where  the  above  assumptions  may  not 
be  appropriate.  Of  particular  concern  is  the  assumption  that  the  failure 
distribution  is  of  the  same  form  at  all  the  stress  levels.  One  reason  for 
this  concern  is  that  different  stress  levels  may  have  different  effects  on 
the  mechanism  which  causes  failure,  and  thus  from  a physical  point  of  view, 
it  may  be  more  realistic  to  allow  for  different  forms  of  the  failure  distri- 
bution at  the  different  stress  levels. 

The  approach  that  we  propose  in  this  paper  requires  neither  distributional 
assumptions  nor  the  specification  of  a time  transformation  function.  Rather, 
our  approach  is  Bayesian,  and  is  prompted  by  what  is  actually  done  in  practice. 
The  Bayesian  point  of  view  allows  us  to  incorporate  some  a priori  information 
which  is  available  in  accelerated  life  tests.  We  would  like  this  paper  to 
be  construed  as  preliminary  and  pragmatic,  and  thus  have  not  attempted  to 
give  a full  theoretical  justification.  Consequently,  we  would  like  to 
invite  mathematical  statisticians  to  resolve  the  statistical  problems  posed. 

1.  Prel iminaries 

As  stated  before,  in  accelerated  life  tests  the  items  are  tests  under 
different  stress  environments.  A stress  environment  may  be  characterized 
by  a single  stress  such  as  voltage  or  temperature,  or  by  multiple  stresses, 
each  of  a different  type.  We  denote  a stress  environment  by  E,  and  the 
set  of  all  E's  by  E.  We  assume  that  the  elements  of  E may  be  ordered 
according  to  the  magnitude  of  their  severity.  Thus,  for  any  two  elements 
E^  and  Ej  belonging  to  E,  E^  > E^  denotes  the  fact  that  E^  is  more 
severe  than  E^.  Let  us  denote  the  k accelerated  stress  environments  by 


Ej.E^ . . . .F^,  and  the  normal  or  use  conditions  stress  environment  by  Eu- 
We  shall  assume  that  the  set  E is  completely  ordered  with  respect  to 
the  relationship  > and  that 


Ej  > ^2  > • • • > E^  > Eu« 


The  basic  problem  in  accelerated  life  testing  is  to  make  inferences 
about  the  life  distribution  of  the  item  under  use  environaent  Eu  using 
failure  data  obtained  under  accelerated  stress  environments  Ej^.  • • • .E^, 
where  k * 1.  In  some  situations,  it  is  possible  to  obtain  a limited  amount 
of  failure  data  under  the  use  conditions  environment  Eu;  however,  in  practice, 
these  situations  do  not  appear  to  be  very  common.  Although  having  failure 
oata  under  the  use  conditions  environment  Eu  has  advantages,  its  absence 
is  in  no  way  detrimental  to  our  approach.  However,  in  order  to  obtain 
results  which  are  useful,  it  is  necessary  that  k be  moderate  to  large, 
and  this  is  what  we  require.  This  requirement  does  not  impose  any  practical 
difficulties  in  many  situations  of  interest,  especially  those  involving 
the  accelerated  life  testing  of  electronic  components  and  in  bioassay  experi- 
ments on  animal  populations. 

In  order  to  introduce  some  notation,  let  us  denote  the  failure  distri- 
bution of  the  items  which  are  tested  under  environment  Ej  by  Fj , where 
Fj(0")  s 0 for  all  values  of  j.  We  assume  that  F^  is  absolutely  con- 
tinuous and  thus  f^(x),  its  probability  density  function,  exists  for 
x e (0,  •).  If  we  denote  1 - Fj(x)  by  F^(x),  then  Xj(x),  the  failure 
rate  of  Fj(x),  is  defined  by 
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fi  (*) 

A.(x)  ■ ^ , F.  (x)  > 0. 

3 Fj(x)  3 


We  find  it  convenient  and  reasonable  to  assume  that  A^(x)  is  continuous 
in  x for  x c [0,  «) . We  use  the  terms  failure  rate  and  hazard  rate 
interchangeably.  Since  Ej  > > • . • > Ek  > E^,  it  is  logical  to  assume 

that 


Aj(x)  i X2(x)  i ...  i Ak(x)  i Au(x) 


for  all  x e [0,  •*) . 

Using  failure  data  obtained  under  Ej,E2,., 

A 

obtain  estimators  A^(x)  of  A^(x),  j - 1,2,., 
0 at  L < » and  all  x e (0,  L], 


(1.1) 


..E^,  we  would  like  to 
.,k,  such  that  for  some 


. st  . st  st  „ 

Aj(x)  i A2(x)  i ...  2.  Ak(x).  (1.2) 

st 

The  notation  X i Y denotes  the  fact  that  X is  stochastically  larger 
than  Y,  i.e., 

P[X£x]  2 P[Y*x]  for  all  values  of  x. 

In  order  to  obtain  estimators  of  A^(x),  j * 1,2, ...,k,  which  satisfy 
Equation  (1.2),  we  shall  use  a Bayesian  approach.  Under  this  approach. 

Condition  (1.1)  is  incorporated  as  a prior  assumption.  Our  approach  is  in 
contrast  with  that  of  Brunk,  Franck,  Hanson  and  Hogg  (l-)66),  who  embody  a 
similar  but  weaker  condition  than  (1.1)  in  their  likelihood  function. 


Specifically,  3runk  et  al.  assume  that  the  distributions  are  ordered,  i.e.. 


Fj(x)  2 F2(x)  2 ...  2 Fk(x)  V x. 


The  above  condition  is  a consequence  of  Condition  (1.1),  but  is  not  equi- 
valent to  it.  Also,  the  approach  of  Brunk  et  al.  is  not  Bayesian. 

We  shall  first  present  a methodology  for  a Bayesian  estimation  of  the 
individual  hazard  rate  yx)  unconstrained  by  (1.1). 


2.  Bayesian  Estimation  of  a Single  Hazard  Rate  Function 

For  estimating  an  individual  hazard  rate  Xj(t),  we  use  the  following 
Bayesian  procedure. 

Let  Nj(t)  be  the  number  of  items  exposed  to  the  environment  Ej  at 
time  t.  Thus,  yO)  is  the  number  of  items  that  are  initially  put  on 
test  in  environment  . 

We  hypothesize  that  the  failure  times  are  governed  by  a time -dependent 
Poisson  process  with  the  probability  of  failure  in  [t,  t*h]  given  by 

yt)A.(t)h  + o(h), 

where  lj(t)  is  the  failure  rate  at  time  t under  environment  E^ . 

For  purposes  of  analysis,  we  divide  the  time  interval  {o,  L]  into 
intervals  of  length  h > o,  where  h is  chosen  to  make  L a multiple  of 
h.  For  convenience,  we  denote  (i-l)h  by  t^,  yy  by  and 

yy  by  A^if  for  i • 1,2, . . . , (L/h) . 


i 


L®t  Xj,i  <*®note  tlle  number  of  failures  in  [t^,  tj+h);  if  there 

are  no  withdrawals,  removals,  censoring,  etc.,  then  N.  ■ N.  . - x. 

J * I J > i j»i 

It  is  helpful  to  clarify  the  above  terms  by  the  following  diagram: 


Number  surviving  * N . . . 


Number  of  failures  * x,  . 

j.i 


oh»t.  lh«t,  2h*t , 

1 * J 


(i-l)h*ti  ih*ti+1 


W hmt (k 


The  joint  distribution  of  xj  f >2' • • • »xj , (L/h)+1»  is  given  by  the 


multivariate  distribution 


(L/h)*l  x.  . / 

ji,!  l:l  A.i1- 


(2.1) 


where  p^^  is  the  probability  of  failure  of  a specified  unit  in  [t^,  t^+p. 


j» (L/h)^l 


(L/h) 

1 " 1 Pji- 

i-1  3>1 


8nd  Xj,(L/h).l“N 


(L/h) 

j»l  ' ih  Xj.i‘ 


Denoting  p^  ^ /E  - by  p^.  Equation  (2.1) 


can  be  expressed 


[Wilks  (1962),  p.  151,  Problem  6.13]  as 


(L/h)  (N  )!  * x.  . * N.  .-x.  . 

A <y;>'  (Vi-'j.i11  (p^)  (1’pi>i)  • 


(2.2) 


The  p*  ,'s  can  be  interpreted  as  the  average  failure  rate  over  the 
3 »i 

interval  [t.,  t.  );  that  is,  p.  . is  the  conditional  probability  that 
1 l+l  3 > 1 

an  item  which  survived  to  time  t^  will  fail  by  time  t^+j.  terms  of 

the  A.  . , we  can  write 
3 »* 


A.  .h  * p.  - ♦ o(h) . 

Jl  1 J 9 A 

liquation  (2.2)  provides  us  with  a starting  point  for  a Bayesian  analysis 


2.1.  Prior  Distributions  for  p.  .. 

J » ~ 

We  shall  confine  the  discussion  in  this  section  to  the  jth  environmental 
condition. 

Suppose  that  the  p*  , i - 1,2, . . . , (L/h) , are  a priori  independent 
3»* 

beta  random  variables  with  marginal  densities 


* r<“i  i*si  * ai  i-1  * si  r1 

f(pM>  - ■ 

J*1  Jl1 

It  is  of  interest  to  note  [cf . Lochner  (1975)]  that  the  above  prior 
density  leads  us  to  a generalized  Diriohlet  density  for  the 

£j  * (Pj,l*pj,2’*-"pj,(L/h)+l)- 

Given  the  N.  .’s  and  the  x.  .'s,  the  posterior  density  of 
3 « 1 Ji1 

* * * 

(pj,l*Pj,2 Pj.a/h)3  is 


- 8 - 


(L/h)  r(a.  .*a  i+N  ) 

n Lt± J » r J > ± 


« rdi : i^i . i-*j  ip-  ^ 


* i^xi  i_i 

\ J 9 1 J 9 1 


(2.3) 


see  DeGroot  (1970)  p.  160. 


3.  A Bayesian  Estimation  of  Ordered  Average  Failure  Rates. 


In  a Bayesian  context.  Condition  (1.1)  leads  us  to  the  requirement  that 


for  every  fixed  value  of  i. 


w w 

>•  i i 2 Pi  i for  j 3 2,3, ... ,k. 

J"1!1  J 9 1 


(3.1) 


w w 

Thus,  our  prior  distributions  on  p.  . . and  p.  . will  have  to 

3 9 1 


chosen  such  that 


ptPi  t , 2 Pi  ^ p[Pi  i 2 P]  *0T  a11  P 2 0. 

J 9 *• 


(3.2) 


One  way  of  achieving  Condition  (3.2)  is  to  assume  that  the  parameters 

* * 

of  the  prior  distributions  of  p.  . . and  p.  . satisfy  the  following  con- 

] "i » 1 J • 1 

ditions  for  every  fixed  value  of  i: 


a.  . s a.  , . 
J.i  J-l.i 


(3.3) 


Bj.i  ^ f°r  j * 2,3, ... ,k. 


For  a proof  of  the  above  statement,  see  Appendix  A. 


r 
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In  order  to  be  assured  that  Condition  (3.1)  is  also  satisfied  with 

respect  to  the  posterior  distributions  of  p*  and  p*  it  is  suf- 

3 "i* 1 iii 

ficient  (see  Appendix  A)  that,  for  every  fixed  value  of  i 


a,  . ♦ x.  . s a.  . . + x.  , 

j.i  3, i J-l,i  j-l.i 


and 


(3.4) 


In  order  to  make  our  Bayesian  analysis  more  practical,  we  will  have 
to  reduce  the  number  of  prior  parameters.  One  way  of  doing  this  is  to  assume 


that 


Oj  ^ ■ a for  all  values  of  i and  j. 


and  that  (3.5) 

B.  . - Bi  for  i - 1,2, , (L/h) . 

J > A J 

Thus,  the  prior  distribution  of  each  p*  , i » 1 ,2, . . . , (L/h) , is  a beta  with 
parameters  a and  B-  This  plus  Conditions  (3.3)  and  (3.4)  lead  us  to  the 
following  remarks. 

The  first  part  of  Condition  (3.4)  will  be  satisfied  whenever  x,  . s x,  . 
for  j » 2,3, ...,k.  That  is,  the  number  of  failures  in  the  interval  (t^,  t^*h) 
under  environment  E^  must  not  be  greater  than  the  number  of  failures  in  the 
same  interval  under  environment  for  all  values  of  j.  Because  Condition 

(3.4)  will  have  to  be  satisfied  for  every  fixed  value  of  i,  i ■ 1,2, . . . , (L/h) , 


- 10 


a reasonable  strategy  is  to  have  N.  . < N.  . . for  j = 2,3,..., k.  That 

J "*i  * 

is,  we  shall  put  more  items  on  test  (initially)  under  the  more  severe  environ- 
ment E.  ,,  than  under  the  environment  E.. 

J-1  J 

Since  N.  . = N.  . . - x.  . . , i = 1,2, , . . , (L/h) , the  second  part 

J > 1 J*1"-1  J > 


of  Condition  (3.4)  can  be  written  as 


0j-l  + 5 0j  + Nj,i+1  for  j * 2'3 k‘ 


(3.6) 


An  interpretation  of  the  above  condition  is  that,  for  every  fixed  value 
of  i,  the  number  surviving  at  the  start  of  the  (i+l)st  interval  plus  the 
prior  parameter  B j ^ for  the  environment  E^  ^ , must  be  smaller  than  the 
corresponding  sum  for  the  environment  E^.  Thus,  whether  the  second  part  of 
Condition  (3.4)  is  satisfied  or  not,  depends  not  only  on  the  number  of  failures 
in  the  it*1  interval  and  the  number  surviving,  but  also  on  the  values  of  the 
pricr  parf-meters  $.  and  Since  the  number  of  failures  in  a particular 

interval  is  a function  of  the  severity  of  the  environment  and  the  number  on 
te-'t,  end  since  ^ has  to  be  less  than  or  equal  to  6^  (see  Conditions 
(3.5)  and  (3.6)),  it  is  reasonable  to  have  0^  ^ < $.  whenever  E^  < E^ 

71. at  is,  the  values  of  the  prior  parameters  0^  and  0^  are  indicative  of 
♦he  relative  severity  of  the  environmental  conditions  E^  ^ and  E^.  Since 
Ej  > E2  > . . . > Ek,  we  will  choose  the  B^'s  such  that  0j  < &2  K * * * < 0k* 
id  the  values  of  the  0^’s  will  be  indicative  of  the  severity  of  the  environ- 
mental conditions  E^,  j * l,2,...,k. 

If  the  prior  parameters  Bj,  j ■ l,2,...,k,  and  the  data  from  the 
accelerated  life  test,  i and  i = 1,2, . . . , (L/h) , are  such  that 


Condition  (3.4)  is  satisfied  for  every  fixed  value  of  i,  then  the  stochastic 

ordering  Condition  (3.1)  will  be  automatically  satisfied  with  respect  to 

* * * 

the  posterior  distribution  of  (p.  .,p.  _,..., p.  ..,.).  If  the  above  is 

J > * Ji^  ji l«/nj 

not  the  case,  then  we  will  have  to  pool  the  adjacent  violators  using  the 
pooling  procedure  described  in  the  next  section. 


4.0.  The  Pooling  of  Adjacent  Violators. 

The  procedure  for  pooling  adjacent  violators  described  here  is  commonly 
used  in  ieotonic  regression)  see  Barlow,  Bartholomew,  Bremner  and  Brunk  (1972). 
The  pooling  is  between  the  violators  of  the  assumed  ordering;  that  is,  whenever 
x.  . . < x.  .,  then  x.  . . and  x.  . are  pooled. 

Consider  the  time  interval  [ (i- 1) h, ih) ; by  Condition  (3.4)  we  require 


xl,i  2 x2,i  * 


iXj-l 


2 x.  . a 


and 


81  + N1  i " X1  i S + N2  i " x2  i 5 5 ®j-l  + Nj-1  i " xj-l  i 


s S.  ♦ N j i " xj  j S ...  5 * Nk  i ” *k  i ‘ 


If  the  reversal  occurs,  that  is,  if  either 


XM.i  * Xj.i 


or  if 


Sj-1  + " Xj-l,i  > Bj  * Nj,i 


- x 


j.i’ 


then  we  pool  the  violators  and  replace  them  as  shown  below. 
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Replace  both  x.  . . and  x.  . by  4<x.  , . ♦ x.  .)  and 

j“i,i  3»*  i j-1,1  j,i 

®j-i  * "j-i.i ' Vi.i  "d  Bj  * Nj,i  • xj,i  bJ' 

* NJ,i  * 

We  now  test  to  see  if  the  new  sequence  is  properly  ordered,  i.e. , 


. i x-  . s ...  2 ^-(x.  . . ♦ x.  .)  ■ i(x.  , . ♦ x, 

L,i  ^»i  2 j-l.i  2 j-l, i j,; 


) * 


and 


Vl.i  * •**  * *lc#i* 


01  * Nl,i  ' xl,i  S s I^j-l  + ej  * Nj-i,i  * Nj,i  " xj-l,i  " xj,iJ 


" + Sj  * + Nj,i  " Xj-l*i  ' Xj.iJ 


s ...  s fi.  + N. 


k + Nk,i  " ^,1' 


If  a reversal  exists  in  either  of  the  new  sequences,  then  we  replace  again 
by  appropriate  averages.  Thus,  if 


^Xj-l.i  * “j.i5  " ^Xj-l.i  + Xj.i}  < Xj+l,i* 


or  if 


2<0j-l  * * Nj-l.i  * Nj,i  ' xj-l,i  ‘ xj,iJ  > 3j*l  * Nj*l.i  " xj*l#i* 


then  we  replace  each  one  of  the  three  by  the  corresponding  average 


frj-M  * *j,i  * 
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L 


(in  the  first  sequence)  and 


**J-»  * *J  * 8j.l  * "j-l.i  * "j.,  * - «,.! 


(in  the  second  sequence). 

We  continue  the  above  procedure  until  all  reversals  in  the  interval 
ih)  ®re  eliminated.  We  use  this  pooling  scheme  for  each  of  the 
(L/h)  time  intervals  to  achieve  the  desired  ordering. 


4*1*  Some  General  Consents  Regarding  the  Pooling  Procedure 
The  following  comments  regarding  the  pooling  of  adjacent  violators  will 
be  helpful. 

1.  An  excessive  amount  of  pooling  occurs  if  the  relationship  specified 

by  Condition  (1.1)  is  incorrect,  or  if  the  environmental  conditions  Ej.Ej Efc 

are  too  similar  to  each  other. 

2.  If  A^(0  ) = 0 for  j « l,2,...,k,  then  it  is  reasonable  to  expect 

that  some  pooling  will  be  necessary  at  the  lower  values  of  i,  i ■ 1,2, . . . , (L/h) . 

For  many  practical  situations,  it  is  reasonable  to  assume  that  the  failure 
distributions  F^,  j * l,2,...,k,u,  have  increasing  failure  rate  (see  Barlow 
and  Pro s chan  (1975)).  If  we  wish  to  make  such  an  assumption,  and  incorporate 
it  into  our  analysis,  then,  for  each  value  of  j,  we  must  have  in  addition  to 
Condition  (3.1), 

} 

Xj,i-1  * Aj,i’ 

or  that 


* st  * 

* ^j,i*  X " 1»2, ....  (L/h) . 


» i 

J 
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The  above  condition  will  further  complicate  our  pooling  procedure;  we 
shall  therefore  not  assume  that  the  are  increasing  failure  rate. 


5.  The  Posterior  Distributions  of  the  p.  . 's. 
lit!— 

It  is  apparent  from  the  discussion  in  Sections  3 and  4 that  the  posterior 
# 

distributions  of  the  p^  ^ 1 s depend  on  the  outcome  of  our  pooling  procedure. 

The  posterior  analysis  is  straightforward  if  no  pooling  is  necessary,  for 

* * * 

then,  the  posterior  distribution  of  (p^  1#p^  2,...,p^  (L/h)-fp  is 


(L/h)  r(a*6.*N.  .)  * a+x.  .-1  * B.+N.  4-x,  4-l 

ill  ,p^] 


for  j • 1,2, ...,k. 

Under  the  assumption  of  a squared  error  loss  function,  a Bayes  estimator 
of  Pj  i ■ 1,2, . . . , (L/h) , j * l,2,...,k,  is  simply  the  posterior  mean; 

that  is 


'j.i  (m+Bj+Nj^)  » 


(see  DeGroot  (1970)  p.  40). 


If  pooling  is  necessary,  then  some  or  all  of  the  x.  ,'s  and  the 

J i* 

(Bj*Nj  j)'s  will  be  replaced  by  their  appropriately  pooled  averages.  In 
any  case,  the  general  expression  for  p.  . will  be  of  the  form  given  above. 

Jl1 


6.  A Model  for  Extrapolations  to  Use  Conditions  Stress. 

Our  analysis  leads  us  to  an  array  of  Bayes  estiaators  of  the  average 
failure  rates  over  intervals  of  length  h,  for  each  environmental  condition. 
Because  of  our  pooling  strategy,  the  Bayes  estiaators  will  be  correctly 
ordered.  That  is,  for  each  fixed  value  of  i. 


pl,i  2 p2,i  2 2 pk,i’ 


(6.1) 


Given  the  p.  .fs,  j ■ l,2,...,k  and  i ■ 1,2,..., (L/h),  our  goal  is 
j # 1 

• • • _ 

to  predict  p .,p  p /M#  the  average  failure  rates  over  the  ti»e 

u,i  u,«  u,  iL/nj 

intervals  [(i-l)h, ih) , i ■ 1,2, . . . , (L/h) , respectively,  under  the  use 
conditions  environment  Eu< 

In  the  absence  of  any  knowledge  about  a relationship  between  the  average 
failure  rates  and  the  values  of  the  various  stresses  which  constitute  an 
environaent,  some  fora  of  an  assuaption  is  essential.  This  is  particularly 
crucial  if  k,  the  nuafcer  of  distinct  environmental  conditions  of  interest, 
is  small.  If,  however,  k is  large,  then  a relationship  between  the  average 
failure  rates  and  the  stresses  can  be  eapirically  obtained;  this  is  what  is 
often  done  in  practice. 

When  k is  small,  we  shall  postulate  the  following  si^>le  but  reasonable 
relationship  between  the  average  failure  rate  estiaates. 

For  some  k unknown  constants  Wq.Wj,. . . ,w^_j,  We  assuaie  *****  *or  each 
value  of  i,  i ■ 1,2, . . . , (L/h) , 


I’m  ‘ “o  * Vk-i.i  * Vk-2.1  • •••  * ’W’l.f 


• a • 


(6.2) 
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The  above  relationship  states  that  the  average  failure  rate  over  a 
particular  tiae  interval  under  the  environaental  condition  is  a weighted 
sun  of  the  average  failure  rates  over  the  saae  tine  interval  under  the  condi- 
tions This  is  reminiscent  of  an  autoregreaeive  prooese 

of  order  k-1  which  has  found  useful  applications  in  forecasting  (see  Box 
and  Jenkins  (1975)). 

In  order  to  sake  the  above  relationship  more  meaningful  we  shall  require 
that  the  environaental  conditions  E^  < j < . . . < E2  < Ej  increase  in 

aagnitude  of  severity  by  the  saae  fixed  amount.  For  example,  if  represents 
a single  stress,  say  a voltage  stress  V^,  then  we  shall  require  that  ^ - Vj  ■ C, 
where  C is  sone  suitable  constant,  j - l,2,...,k. 

Since  Equation  (6.2)  holds  for  i ■ 1,2, . . . , (L/h) , the  least  squares 
estimators  of  w0,wi,****,l'k-l  can  be  °'>taine^  in  a routine  manneT.  These 

A A A 

estimators  are  denoted  by 

If  the  increase  in  the  magnitude  of  the  severity  of  the  environment  from 
Eu  to  E^  is  the  same  as  that  from  E^  to  E^j,  j • 2,3, ...,k,  then  an 
estimator  of  the  average  failure  rate  under  Eu  is  given  by 


A • 

P 


u,i 


* Vk.i  * 


♦ w, 


'k-l,ip2,i 


(6.3) 


(L/h)  .. 

for  i » 1,2,..., (L/h),  and  Pu>(L/hj+1  * 1 “ ^ Pu,i* 

If  the  increase  in  magnitude  of  the  severity  of  the  environment  from  Eu 
to  E^  is  two  times  that  from  Ej  to  E^j,  j ■ 2,3,...,k,  then  we  iterate 


I 
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> 

I 


upon  Equation  (6.3)  one  more  tiae  to  obtain  the  desired  result.  Thus,  in 
principle,  we  can  iterate  upon  Equation  (6.3)  as  many  tines  as  is  necessary, 
depending  upon  the  separation  between  Ey  and  E^. 

* 

By  the  definition  of  Pu  the  probability  of  an  itea  surviving  to  tiae 
t (assuaed  a aultiple  of  h)  under  environaent  Ey  is 


_ * t*/h  * 

-uCt  ) - £ d-Pu>l). 


Thus,  we  have  the  following  as  an  estiaator  of 


(t  ): 


« • t*/h  » * 

Fu(t  ) - £ 


where  the  p . are  given  by  (6.3). 

U«  1 

The  properties  of  this  estiaator  have  not  been  studied.  The  estiaation 
procedure  has  been  proposed  as  a practical  method  for  providing  answers  in 
the  difficult  area  of  accelerated  testing. 


APPENDIX  A. 


In  this  appendix  we  determine  sufficient  conditions  under  which  two  beta 
random  variables  are  stochastically  ordered.  The  results  may  already  be  known; 
however,  they  are  included  for  completeness. 

The  notation  and  terminology  are  taken  from  Barlow  and  PToschan  (1975). 

xa_1e'x 

We  first  note  that  a gamma  density  f(x,  o)  « — f^)-  is  TP2  in 

* ya‘l  _x 

(x,  a).  Therefore,  F(x,  a)  » / jrr^r  e dx  is  TP2  in  (x,  o)  since 

x 1 J 

F(x,  a)  « / f(y,  a)H(x,  y)dy,  where  H(x,  y)  » 1 for  y * x and  H(x,  y)  ■ 0 


y < x.  (Note  that 


IhUj,  yx)  H(x1,  y2) 


0 for  x.  < x-,  y.  < y0.  Thus, 


Imxj,  yx)  H(x2,  y2)| 

H(x,  y)  is  TP2,  and  so  the  composition  of  f,  H is  TP2>) 
Since  F(x,  o)  is  TP2,  then  for  xJ  < x2,  < a^, 

lf("l-  V F(Xl’“2,L„  ~ „ I * 


|F(x2,  Oj)  F(x2,  a2)| 


i.  0.  Choose  Xj  * 0.  Then 


|F(x.  Oj)  F(Xf  <*2) 


for  0 < x < «»,  Oj  < Oj.  Thus,  F(x,  a2)  2 F(x,  ap  for  and  x > 0. 


We  have  shown: 


1.  Lemma.  Let  X be  a gamma  r.v.  with  shape  parameter  a.  Then  Xq 


is  increasing  stochastically  in  a. 


2.  Proposition.  Let 


Y«t  have  density  ^ 


V1 

hr- t 1 * 1>2>  and  • 


Y be  independent.  Then  X 


has  beta  density 
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r(<Va2}  “r1  a2_1 

f(Oj,a2,x)  - 7fr“)f(a2)'  * 0-x)  for  OSxU.  (Hogg  and  Craig,  1970, 

p.  134). 


3.  Tlworen.  Let  X have  beta  density  f(a.,a,,x).  Then  X 

1*  2 1 1 °i*°2 

is  increasing  stochastically  in  and  decreasing  stochastically  in  . 

Ya  Yn 

Proof.  Write  X - 


- 1 v _ 1 • 
a « y «,  v < \ -i  * v I“v — * where  a.  < a-. 

V 2 ‘aj  la2  “l*°2  Y0l  * V 22 

For  fixed  By  unconditionin*  X^  ^ , for 


®2  < a2* 


st 

By  sinilar  reasoning,  wc  can  show  X s x , for  a < o'. 

V°2  al’a  2 1 1 


REFERENCES 


[1]  BARLOW,  R.  E.,  D.  J.  BARTHOLOl'EW,  J.  M.  3REMNER  and  H.  D.  BRUNK  (1972). 

Statistical  Inference  Under  Order  Restrictions.  Wiley,  New  York. 

[2]  BARLOW,  R.  E.  and  F.  PROSCHAN  (1975).  Statistical  Theory  of  Reliability 

and  Life  Testing.  Holt,  Rinehart,  and  Winston,  New  York. 

[3]  BOX,  G.  E.  P.  and  G.  M.  JENKINS  (1976).  Time  Series  Analysis,  Forecasting 

and  Control.  Holden-Day,  San  Francisco. 

[4]  BRUNK,  H.  D.,  W.  E.  FRANCK,  D.  L.  HANSON  and  R.  V.  HOGG  (1966).  Maximum 

likelihood  estimation  of  the  distributions  of  two  stochastically 
ordered  random  variables.  J.  Aroer.  Statist.  Assoc.  6 1067-1080. 

[5]  DeGROOT,  M.  H.  (1970).  Optimal  Statistical  Decision.  McGraw-Hill, 

New  York. 

[6]  LQCHNER,  R.  H.  (1975).  A generalized  Dirichlet  distribution  in  Bayesian 

life  testing.  J.  Roy.  Statist.  Soc.  B 37  103-113. 

[7]  MANN,  N.  R.,  R.  E.  SCHAFER,  and  N.  D.  SINGPURWALLA  (1974).  Methods 

for  Statistical  Analysis  of  Reliability  and  Life  Data.  Wiley, 

New  YorlT 

[8]  PATHAK,  P.  K.,  W.  J.  ZIMMER,  and  R.  E.  WILLIAMS  (1977).  Estimates 

of  parameters  in  acceleration  models.  Technical  Report  No.  334. 

The  University  of  New  Mexico,  Albuquerque,  New  Mexico. 

[9]  SHAKED,  M. , IV.  J.  ZIMMER,  and  C.  A.  BALL  (1977).  A nonparametric  approach 

to  accelerated  life  testing.  Technical  Report  No.  340.  The  University 
of  New  Mexico,  Albuquerque,  New  Mexico. 

[10]  WILKS,  S.  S.  (1963).  Mathematical  Statistics.  Wiley,  New  York. 


- 20  - 


